Anharmonic transitions in nearly dry L-cysteine I 
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Two special dynamical transitions of universal character have been recently observed in macro- 
molecules at Td ~ 180 — 220 K and T* ~ 100 K. Despite their relevance, a complete understanding of 
the nature of these transitions and their consequences for the bio-activity of the macromolecule is still 
lacking. Our results and analysis concerning the temperature dependence of structural, vibrational 
and thermodynamical properties of the orthorhombic polymorph of the amino acid L-cysteine (at a 
hydration level of 3.5%) indicated that the two referred temperatures define the triggering of very 
simple and specific events that govern all the biochemical interactions of the biomolecule: activation 
of rigid rotors (T < T* ), phonon-phonon interactions with phonons of water dimer (T* < T < Td), 
and water rotational barriers surpassing (T > Td). 
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Recently two special dynamical transitions, apparently 
of universal character, were observed in macromolecules 
as lysozyme, myoglobin, bacteriorhopsin, DNA, and 
RNA (see, e.g., refs.p], Q). With cooling, the first tran- 
sition occurred at Tjj « 180 — 220 K at hydration level 
h > 18%. It is related to a transition from anharmonic 
to harmonic regime [2]. The most remarkable character 
of this transition is its direct correlation to the biological 
activity p], Q. On the basis of quasi-elastic neutron scat- 
tering (QENS) measurements on lysozyme, Chen et al.[3| 
interpreted this transition as a fragile-to-strong dynamic 
crossover where the structured water makes a transition 
from a high-density to a low-density state. This scenario 
was based on studies that advocated a second critical 
point of water (liquid-liquid critical point) in confined 
supercooled water [4]. However, a recent high resolution 
QENS experiment performed by Dosteret al.[5j in fully 
deuterated C-phycocyanin protein, showed no evidence 
of such fragile-to-strong character in the dynamical tran- 
sition at Td. As pointed out by these authors, their ex- 
perimental findings are consistent with a glass transition 
of the hydration shell. However, the specific nature of 
the glassy state remains unclear. Moreover, the specific 
role of protein and water is still an issue. As an example, 
Frauenfelder et al.[l] proposed that the protein dynamics 
is slaved to the bulk solvent fluctuations. On the other 
hand, the dynamics of tRNA compared to other studies 
on DNA and proteins by Khodadadiet al.[6| contradicts 
the slave picture. Their results indicated the mutual in- 
fluence of macromolecule and hydration water must be 
considered. The other important transition observed in 
macromolecules is a hydration independent sub-regime 
of harmonicity onset at T* ~ 100 K[2]. The microscopic 
nature of this transition is still an object of debate. The 
discussed possibilities include quantum effects near a zero 
point vibration[7j, and methyl group rotation associated 
with anharmonicity[2[ Q. Vibrational modes play a very 
important role in the conformational arrangement of dif- 
ferent functional states and also in the energetic balance 
of the chemical bonds [9]. Theoretical studies on the an- 
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FIG. 1: Experimental (open circles) and calculated (black 
solid line) diffraction pattern of orthorhombic L-cysteine I 
structure. The vertical bars and red solid line at the bottom 
represent the Bragg reflections and the difference between ob- 
served and calculated data, respectively. The insert shows the 
unit cell of L-cysteine I observed along the ac (panel I) and 
be (panel II) planes. 



harmonic effects on phonons show that their frequency 
and linewidth analyzed by the inelastic light scattering 
(Raman effect) could vary due to the anharmonic poten- 
tial shape of the simple oscillator and also due to phonon- 
phonon interactions | lOl [Tlj . 

The present study is focused on the temperature de- 
pendence of structural, Raman-active vibrational modes 
and thermodynamical properties of the orthorhombic 
polymorph of the amino acid L-cysteine, C3H7NO2S 
+NH 3 - CH(CH 2 SH) - COO' . This amino acid pos- 
sesses a very simple chemical structure and high bio- 
logical relevance. The thiol or sulfhydryl group in the 
residues of L-cysteine is the most chemically reactive site 
in proteins under physiological conditions [12|. Therefore, 
the knowledge of its dynamics is very important. 

Figure 1 shows the X-ray powder diffraction pattern of 
the studied sample. The Rietveld refinement indicated 
that the predominant phase is orthorhombic {P2\l\l\ 
space group) with unit cell parameters a = 8.09758(8), 
b = 12.1495(1), and c = 5.4142(1) A. The unit cell (in- 
sert of Fig. 1) has 4 molecules (2 = 4) of L-cysteine. Our 
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FIG. 2: Raman spectra at 30 K (black line) and 300 K (red 
line). Those bands marked with asterisk (*) had their detailed 
temperature behavior. 



findings are consistent with others previously reported in 
the literature [HQ. 

Figure 2 shows the Raman spectra of L-cysteine I at 
30 and 300 K. A clear energy threshold at 500 cm -1 (60 
meV) could be observed. All bands above this energy 
kept their shape with merely a softening due to tem- 
perature variation. The exception was the bands at 815 
and 830 cm -1 where an intensity transfer was observed 
between them. A completely diverse situation occurred 
below 500 cm -1 where some bands broadened at high 
temperatures. Ab initio vibrational analysis calculations 
were performed to interpret the Raman spectra. Since 
the relevance of the intermolecular interactions (cys-cys, 
H-bond, etc) is an important issue, two methodologies 
based on the Density Functional Theory (DFT)[15] were 
compared: free molecule calculations performed using 
the B3PW91 hybrid functional[16l as implemented on 
the on Gaussian 03 software suite [13 and unit cell with 
periodic boundary condition calculations using a fixed 
unit cell volume and the BLYP functional [l8| augmented 
with corrections for the proper description of dispersion 
interactions [19] as implemented in the CPMD code[20|. 
Only those values in accordance with the experimental 
finding within 5% variation were considered for the as- 
signment. It is important to stress that the free molecule 
calculated frequencies that meet this criterion were only 
those modes related to carbon backbone torsion. Similar 
results were obtained by Pawlucojk et al.|l4[. The cor- 
responding bands in the crystal calculation showed bet- 
ter agreement with experimental values with a larger set 
of concordant bands. As a matter of fact, most crystal 
modes presented contributions of almost all atoms in the 
unit cell with their corresponding eigenvectors showing 
a more complex dynamics than the free molecule case. 
This fact is a manifestation of the role of the intra and 
inter-molecular interactions in the vibrational modes de- 
scription, both responsible for the strong anharmonicity 
of this system. 

We considered a detailed temperature dependence 
study for those modes where hydrogen bonds participate 
(signed with asterisks on Fig. 2). These particular vi- 





„ O G1 






b) 

" O G2 






□ G3 




1007 cm 1 








FIG. 3: % a oj/oj (black scale) scaled to (% A V/V ) (blue scale) 
as a function of temperature for Gl (open circles,a)), G2 
( open diamond, b)), and G3 (open square, b)) phonons. 
(%ACo> /oJcorr) is shown in the right side (c) and d)). The solid 
lines show the fitting to the corresponding polynomial func- 
tions (quadratic, c); and linear, d)). The inserts represent the 
unit cell motion evolved in each phonon group. 



brations are candidates to mediate the water-molecule 
interactions. At this point it is important to distinguish 
the two kinds of anharmonic contributions to the temper- 
ature dependence of the phonon frequency. The first one 
arises from anharmonic deviations of the single harmonic 
oscillator potential shape. The second one is related 
to phonon-phonon coupling [llj . The isothermal mode 
Griineisen parameter is a parameter that quantifies the 
phonon anharmonicity of the first type. A good estima- 
tion for could be obtained scaling the phonon frequency 
(Auj/uj ) to unit cell (AV/V) temperature variation. 

This scaling was performed on the phonons indicated 
by asterisks. We were able to distinguish three kinds of 
characteristic behaviors. Figure 3a) and b) display the 
representative behavior for each member of these groups. 
The right and left scales show the phonon frequency vari- 
ation and the negative of the unit cell volume variations, 
respectively. The first group (Gl) comprises the 455, 874, 
and 945 cm -1 phonon modes. Figure 3a) (open circles) 
shows the representative temperature dependence for the 
455 cm -1 vibration. With heating, an overall softening 
for all members of this group following the lattice expan- 
sion was observed. However, the softening rate displayed 
a smooth increase between ~ 100 and 250 K compared 
to AV/V. The second group (G2) comprised the 538, 755, 
and 771 cm -1 phonons. The observed behavior was the 
opposite of the one observed following lattice expansion, 
as can be seen in Fig. 3b) (open diamond) for the 771 
cm -1 mode. The third group (G3) has a single member 
(1007 cm -1 phonon) and presented a very peculiar be- 
havior (open square symbol in Fig. 3b). Up to 70 K, 
its frequency was almost constant. Thereafter, a sudden 
hardening was observed between 70 and 150 K followed 
by a softening. In order to carry out a detailed analy- 
sis of this contribution we subtracted the scaled behavior 
from the phonon temperature dependence. The data is 
shown in Fig. 3c) and d). Those modes of Gl group 
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presented a characteristic quadratic behavior (Fig. 3c). 
The data for G2 and G3 group phonons displayed a clear 
linear dependence. Above ~ 200 K a tendency to change 
was observed for G3, which presented smaller slope than 
the G2 one. It is important to notice that recent reports 
on literature claimed for the occurrence of signatures of 
subtle phase transitions at ~ 70 K in L-cysteine by spe- 
cific heat, Raman spectroscopy, X-ray diffraction, and 
neutron scattering [2 ll - l23j .This kind of behavior can be 
understood by considering phonon-phonon interactions. 
Balkanski, Wallis, and Haro[10| developed a detailed the- 
ory of the higher order phonon anharmonic decay in or- 
der to explain the temperature dependence of the fre- 
quency and linewidth of the optical phonons probed by 
Raman scattering. Since the Debye temperature (0d) of 
L-cysteine is very low (< 5 K as estimated by our specific 
heat data in Fig. 4), we could use the high temperature 
limit of the Balkanski theory. Using this approximation, 
the decay of one phonon with frequency in two and three 
phonons will provide specific temperature contributions 
given by 

Auj/uj oc T, for a three phonon decay process, and 

Auj/oo oc T 2 , for a four phonons decay process. 

Thus, in order to explain the observed temperature 
behavior of the three afore-mentioned groups one needs 
to find phonons that fit the sum rule and would par- 
ticipate in the decaying process. In the present con- 
text the water vibrational bands are natural candi- 
dates. In fact, our vibrational analysis calculations of 
water dimer showed that water dimer presented a set 
of low frequency intermolecular vibrations at 93 cm -1 
(donor twisting/acceptor rocking); 161 cm -1 (donor 
rocking/acceptor wagging); 195 cm -1 (acceptor rocking); 
and 251 cm -1 (0-0 stretching) which could participate in 
the decay process. For Gl phonons, we can establish the 
following correlation 

455 ^93 + 161 + 196 cm" 1 

874 455 + 251 + 161 cm" 1 

945 -+ 455 + 2 x 251 cm -1 , which is correct within 2% 
precision. Each phonon decays into a combination of wa- 
ter and/or Gl-member phonon. Likewise, the correlation 

538 ^ 455 + 93 cm" 1 

755 ^ 538 + 195 cm" 1 

771 -+ 538 + 251 cm -1 , which is correct within 3% 
precision could be established for G2 phonons. Finally, 
for G3 (1007 cm" 1 ), we have 1007 -+ 755 + 251 cm" 1 , 
correct within 2% precision. 

The eigenvectors of the vibrational modes for each 
group also presented characteristic patterns. The vi- 
brational modes of Gl group were characterized by a 
stretching motion of the unit cell along the three crys- 
tallographic axes (see inset of Fig. 3c). The atomic mo- 
tions in the G2 and G3 groups resulted in the stretching 
of the cell along two perpendicular axes (see inset of Fig. 
3d)). For G2, the movement was along b and c direc- 
tions. In the G3 mode the movement was confined in 
the ab plane. This symmetry may be relevant to the un- 
derstanding of the interactions between L-cysteine and 



water modes. In fact, our results concerning the water 
dimer vibrational modes indicated that the dipole deriva- 
tive unit vector lied along the donor oxygen and acceptor 
hydrogen direction (main component along the 0-0 axis) 
for 93 and 161 cm -1 modes, being restricted to the per- 
pendicular directions for the other modes. Changes in 
the relative orientation of the hydrogen atoms of the wa- 
ter dimers might have a strong influence on the coupling 
to the L-cysteine modes. The rearrangement of the hy- 
drogen bonding pattern due to thermally activated tun- 
neling pathways among configurations of water dimers 
with exchange of acceptor and donor hydrogens [2J is an 
option that should be considered. Our results could be 
consistently explained once considering the presence of 
one water dimer in the L-cysteine unit cell. Since the 
G2 phonons do not change their behavior above 100 K 
(Fig. 3d) we argue that the dimers have their 0-0 axis 
along the a-direction. At 200 K, there is enough ther- 
mal energy to overcome the three water dimer tunneling 
pathways (bifurcation, interchange, and acceptor switch- 
ing, see Keutsche al. in ref. [24]. This might explain 
the slope change observed on the G3 phonon frequency 
behavior (Fig. 3d). 
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FIG. 4: a)C p showing two tiny transitions at T* ~ 80 K and 
Td ~ 230 K. The solid lines are simulations as described in 
the text. 

The specific heat (C p ) results are shown in Fig. 4. Fig. 
4a) shows the data between 50 and 300 K. Very sub- 
tle transitions are seen at 220 K and 80 K as broad 
and tiny bumps near these temperatures. Calorimet- 
ric measurements on orthorhombic L-cysteine in the lit- 
erature reported only a subtle transition at T* = 70 
K[2l[. Incoherent inelastic neutron scattering on mon- 
oclinic L-cysteine [23| revealed an anharmonic transition 
at Td = 150 K. Thus, it is reasonable to associate the 
two transitions observed for L-cysteine to the Td and 
T* ones, which is commonly seen in proteins and other 
macromolecules . 

According to the literature, the low temperature tran- 
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sition at T* is related to methyl group rotations in 
proteins Since the specific heat of anisotropic rigid 
rotors is well known (25j. this hypothesis could be easily 
checked. Thus, considering additional linear (electronic) 
and cubic (phononic) contributions, the specific heat is 
written as 

C p = 7 T + f3T 3 + C rotor (T, I xy , I z ) 

where C ro tor(T, I xy , I z is given by expressions (3) and 
(4) of ref.[25], with I xy and I z the inertia momenta in the 
xy— plane and z— axis, respectively. The simulated curve 
is shown as a solid line in Fig. 4a). The set of parameters 
that furnished the best simulated curve for T < 150 K 
were found to be 7 = 0.205 J/mole K 2 ;f3 = 1.15 x 10~ 6 
J/mole K 4 ; and I xy /I z = 0.72 . The high 7- value is 
consistent and expected for an isolating material. The 
Debye temperature estimated from j3 value was Qd < 
5K. The ratio I xy /I z found is consistent with a prolate 
symmetric top rotor (CH2 rotation). 

It is a remarkable fact the observation of the Tjj tran- 
sition at very low water content (3.5%) as described 
herein. This water concentration is equivalent to 1 water 
molecule per unit cell in average. This finding is incon- 
sistent with the fragile-to-strong picture of Chen et al[3| 
, since it is not possible to treat the aqueous content as 
very low density aqueous phase. In a similar way, the 
idea of a glassy transition of the hydration shell 0, H| 
is not supported. Thus, another explanation concern- 



ing the origin of the dynamical transition should be con- 
sidered. By subtracting the above-mentioned simulated 
curve from the raw data, as it can be seen in Fig. 4b), 
we obtained an estimate of the heat capacity evolved in 
the transition To. The broad transition resembles the 
shape of a Schottky anomaly due to thermal popula- 
tion in some specific energy levels. One possible inter- 
pretation is to consider the anomaly as originated from 
thermal population of tunneling splitting energy levels in 
configurations of water dimers to be consistent with the 
analysis of Raman data in Fig. 3d). The best choice 
was obtained using the three-level scheme with energies 
Ei = 250K (174 cm" 1 ), E 2 = 600 K (420 cm" 1 ), and 
E 3 = 900 K (629 cm -1 ). The agreement between the 
experimental data and simulation was reasonable consid- 
ering the data dispersion. Keutsch et al.[24] calculations 
(gaseous phase) indicated that the three distinct water 
dimer low barrier tunneling pathways acceptor switching, 
interchange tunneling, and bifurcation tunneling have en- 
ergies of ~ 157, 207, and 394 cm -1 , respectively. The ex- 
tra constraints imposed by the environmental hydrogen 
bonds will increase the tunneling splitting. Therefore, 
the energies obtained in the simulation compared to the 
gaseous phase dimer are expected to be greater. 

This work was supported by Fapesp/Brazil, 
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